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Abstract 

A novel iterative algorithm for the efficient computation of the intersection areas of an arbitrary number 
of circles is presented. The algorithm, based on a trellis-structure, hinges on two geometric results which 
allow the existence-check and the computation of the area of the intersection regions generated by more than 
three circles by simple algebraic manipulations of the intersection areas of a smaller number of circles. The 
presented algorithm is a powerful tool for the performance analysis of wireless networks, and finds many 
applications, ranging from sensor to cellular networks. As an example of practical application, an insightful 
study of the uplink outage probability of in a wireless network with cooperative access points as a function of 
the transmission power and access point density is presented. 

I. Introduction and Problem Statement 

The computation of the intersection area of many circles is a challenging problem. While the 
intersection of two circles is straightforward, even three circles admit several configurations, each 
resulting in a different expression for the intersection area. Given the centers and the radii of the 
circles, the automatic discrimination among the various cases requires involved condition testing. If we 
consider cases with several circles the problem may appear unsolvable, as the close-form expressions 
for the intersection areas become more and more involved and depend on the specific configuration 
among a huge number of possibilities. 

Despite the wide range of applications of this geometric problem, a systematic approach is 
still lacking. It has been addressed in the literature for three circles, but only for some specific 
configurations, in The intersections among d-dimensional balls are used to find their union in [|2l, 
in a way analogous to the last step of our proposed algorithm, whereas [|3l analyzes the special case 
of equal circular disks, showing that their intersection can be derived from intersection areas among 
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fewer circles. However, no algorithmic solution is proposed to exploit the result in an organized and 
exhaustive way. 

In this paper, an algorithm that efficiently computes the intersections of an arbitrary number of circles 
is presented. The algorithm works in an iterative fashion and is based on a trellis structure. At each 
iteration, the existence of any intersection is checked based on the areas computed in the previous steps, 
thus highly reducing the computational load. Moreover, only the first three steps involve geometric 
considerations, whereas, when the number of circles is higher than three, all the areas can be found 
via simple algebraic calculations. The presented algorithm allows to efficiently solve configurations 
with many tens of circles, without any assumption on the centers and radii of the considered circles. 

The technical contributions of this paper are as follows: 

• we derive two theorems, that provide an easy way to check the existence and calculate the area 
of the intersection region of more than three circles, once the existence and the area of the 
intersections involving a smaller number of circles are known; 

• we present a trellis-based iterative algorithm that allows an easy computation of the wanted areas 
even for configurations with a large number of circles. 

In the following, the applications of the presented tool in wireless networks' performance analysis 
are discussed, and then the geometric problem is formally defined. 

A. Applications to Wireless Networks 

Many frameworks for the evaluation, analysis and simulation of wireless communications are 
based on a signal propagation model in which the attenuation incurred by a transmitted signal is 
a monotonically decreasing function of the distance from its source. Thus, the performance of a 
receiver is a function of its distance from the source, and this leads to a characterization of wireless 
networks based on the concept of coverage range. 

The coverage range of a transmission can be defined by assigning a threshold bit error rate (BER), 
packet error rate (PER) or signal-to-noise-ratio (SNR) which determines an admissible region of 
received power. The coverage range then is the maximum distance between two nodes which guarantees 
the received power to lie within the admissible region. The resulting coverage area of a transmitter 
(receiver), is a circle centered on the receiver (transmitter) and with radius equal to the coverage 
range. Note that different transmission power levels, packet encoding rate and, in general, transmission 
parameters can be represented as multiple circles centered on the same node. 



A node placed in a point of the plane covered by multiple coverage areas can communicate with all 
the nodes associated with those coverage areas. The computation of the area of those regions enables 
a wide range of considerations in many scenarios of interest. 

In cellular networks, circular coverage areas may be used to design base station positioning in order 
to guarantee connectivity [4J. Recently, considerable attention has been devoted to the study of relaying 
strategies in multihop cellular networks to improve capacity, coverage range and Quality-of-Service 
fairness (e.g., see [l5l-[I3). Circular coverage areas of base stations and relays Q can be used to build a 
simple connectivity model aimed at the calculation of the overall capacity of the cell [8J. The ability to 
calculate the area of the various intersections of the coverage areas granted by the proposed algorithm 
may be used to compute the probability that a mobile falls within coverage of a certain set of base 
stations/relays. 

The computation of the intersection areas may also be used to model connectivity in many other 
infrastructured network scenarios. For instance, in heterogeneous networks, the areas covered by 
different network infrastructures (GSM, UMTS, local area networks, and so on) may intersect. Thus, 
the areas covered by multiple technologies may be used in order to allocate users and compute the 
average performance. This problem has been recently investigated in [191 for downlink ^-tier cellular 
networks. 

In non-infrastructured ad hoc networks, circles have been traditionally used to characterize channel 
sensing and data packet decoding. Again, given a topology, the area of the regions in which a new 
transmitter detects/decodes signals from the various sources can be computed using the proposed 
algorithm]^ 

In sensor networks, localization relies on the reception of beacons sent by nodes whose positions 
are known. The accuracy achieved by the localization algorithm depends on the number of beacon 
sources that the node can hear. This requires the computation of the probability that a node falls within 
an area covered by a certain number of circles. Furthermore, intersections of multiple circles are also 
found when addressing the problem of preserving complete sensing coverage of a certain area and 
connectivity ifTOll. [fTTll. 

Another important example in which intersection areas are a fundamental aspect of the performance 

'Multiple circles associated with each base station/relay can be used to account for coverage shrinking as the number of mobiles 
increases. 

^For instance, the algorithm can be directly applied to computing the probability that a node with uniform spatial distribution falls 
within a region connecting other nodes, or becomes a hidden or exposed terminal. 



analysis is routing [12]. The intersection of the circles may represent the area in which a user can 
provide connectivity to some nodes of the network (corresponding to the centers of the various circles). 
When considering geographic packet forwarding [fOI . [[T?]. intersection areas may be helpful to derive 
the distribution of the advancement and the success probability of the communication. 

In this paper, as an example of application of the presented tool, we study the uplink outage 
probability in a wireless network with cooperative access points as a function of the transmission 
power and access point density. An analogous scenario can be found in cellular networks, where 
recent work showed that cooperation among Base Stations may offer considerable performance gain. 
Multi-cell processing (MCP) has been proven to grant higher throughput and achievable data rate [fTSl . 
lfT6l . depending on the topology as well as on the robustness of the backhaul links. Capacity may 
be also increased, as was shown in [[TtI . when cooperation is aimed at cancelling interference. The 
use of relays, together with cooperating Base Stations, has been also considered in [18]. It is clear 
that the performance of MCP depends on how many access points are able to receive and decode the 
transmission from a given source, which can be statistically determined by finding the intersections 
of the coverage areas of Base Stations and relays. 

In the following section, we state in detail the addressed geometric problem and the contribution of 
the paper. 

B. Problem Statement and Contribution 

Consider a set C = {71, 72, . . . , 7Ar^} of circles, whose centers and radii are known. The circles 
in C may partially overlap. We denote with 

-{2^{nl.,„}' ^1' • • • ' ^-^{1, . . . , iVj, z,^z„, forjyn} (1) 

the set of all the possible intersection regions generated by n circles, where x|"^^ .^|= niG{n in} ^« 
the set of the points that belong to all circles in {7ji, 7i2, • • • , 7j„}^C. The set I^^^ contains the circles 
in C. We also define the notation X^") (ii, • • • ,iN^-n) to denote the intersection of n circles out of the Nf. 
in C where circles ii, . . . , iN,-n are not considered, i.e., X^^^n, . . . , iN,-^n)= r\ie{i,...,N,}\{iu-,iN,-„} 
which is hence equivalent to 2:{"^...,;v4\{n,...,i^_„}- 

However, these intersections are not disjoint regions of the plane. See for instance Fig. [T] where a 
configuration with three circles is depicted. In the figure, xj'^^2 3}=-^!' 2r|^''2|=^iU^4, x|^''3|=^iU^2 
and X^^^^^y=AiUA3. 



We call the regions Ai in the figure exclusive intersection regions, as they correspond to the 
intersection of a certain subset of circles, excluding the regions covered by the other circles in C. 
We denote these regions as £^^'\ii, . . . ,iNc-n), where 

ee{ii,.--,«JVc-n} 

For instance, £^^^^(2,3) is the region of the plane covered by 71 and that does not overlap with any 
other circle of C {A5 in Fig. [B, and £^'^^^(3)=^4 is the intersection of 71 and 72, excluding the area 
covered by 73. We define the set as the set of all the exclusive intersection regions generated by 
n circles. Let us define the set £= |Jn=i Af^ ^^^^ ■ Then, the elements of £ are disjoint regions that 
tessellate the overall region covered by C. We also define the measure (p, where (p^A) corresponds to 
the area of the region A. 

In this paper, we address the problem of computing the measure of the regions £^'^\ n=l, . . . , Nc, 
given the centers and the radii of the circles in C. This appears to be an extremely complex geometrical 
problem. In fact, while the area covered by the intersection of two circles has a simple measure, even 
the intersection area of three circles has a rather involved form, that depends on the mutual positions 
of the circles [|T]. When more than three circles are considered, the number of configurations grows 
larger, and the complexity of the geometric conditions and the associated expression of the intersection 
area become difficult to handle. 

The rest of the paper is organized as follows. Section |ll] presents the geometric results that are the 
foundation of the iterative algorithm. In Section |lll] we describe in detail the structure of the proposed 
algorithm. In Section IIII-CI we show how the wanted areas can be computed. Finally, in Section |V] we 
present and discuss the aforementioned network design problem. 

II. Geometric results 

In this section, we present the geometric results that represent the core of the proposed algorithm. A 
key observation is the following: a necessary (but in general not sufficient) condition for the existence 
of the intersection of n circles is the existence of the intersection of all subsets of n — 1 circles]^ This 
consideration may be very useful, since the calculation of several areas among circles which are not all 
intersecting can be avoided. Nonetheless, a stronger result can be stated if the number of considered 
circles is greater than 3. 



^We say that an intersection exists if it is non-empty. 



Theorem 1 

Consider a subset of n<Nc circles 5={7i^ , • • • , 7i„ and the associated intersection region xj"^'' 
With a slight abuse of notation, we refer to the considered subset of circles when denoting the intersection 
regions, and we write Z{^~'li^y\{,^_,^y=^^''-'Kei, • • • , ee). 

Then, V{ii, . . . , in}, ij'^iu, forj^u, if n>A the following holds: if 



{Z^^-'\k))>0 VA;G{zi,...,^4, (3) 



then also ■ x) > 0. 



In other words, the existence of the n intersections among all the subsets of n — 1 circles in S, besides 
being an obvious necessary condition for the intersection among the n circles to exist, is also sufficient 
when n > 4. 

Proof: Consider a set of n circles {7,^, . . . , 7j^}. By hypothesis (j){X^'^~^\k))>0, WkE{ii, . . . , in}- 
Fix the index k=k and define the set 5fc={ii, . . . As 0(X("-i)(I)) > 0, then I^'^-^^^j^d). 

Call A the polygoj^ delimiting X("~i) (/c) and whose sides are m arcs of circumference aj,J=l, . . . ,m, 
with l<m<2{n—2). Note that more than one arc aj may belong to the same circle 7^0 We denote 
with a.{i) the set of arcs belonging to the circle 7^. 
We have to consider three cases: 

1) 37,g5^: aW=0; 

2) V7ie5^|a(?)| = l; 

3) V7.g5^ a(O^0 and 3j : |a(j)|>l 

In the first case 7^ fully contains the whole intersection but its circumference does not hit it. Thus 
j(n-i)(^^)=j(n-2)(^^ z)n7i=X("-2)(^, i) . Slncc by assumption X("-i)(A:)7^0, yke{ii, in}, then also 
X("-i)(i)=X("-2)^^^^) p^__^0_ Therefore, X("-i)(fc)n7^^0, and thus x{"^^ ^.^^^0 and 0(x|"^^ ^.^^)>0. 

In the second case, the polygon A is delimited by exactly n— 1 arcs of circumference, each belonging 
to one of the n—1 circles in S-^. This situation is depicted in Figured Consider two non-consecutive 
arcs of A, namely a-r, at, and assume, without any loss of generality, that they belong to the circles 
7i^ and 7^2, with k^ii, Z2. By assumption, as the intersection of any combination of n—1 circles exists, 
then 7^ must contain at least one point PGX*^"^^^(ii, k) and one point geX{«-2)(^2, k). Since a circle 
is a convex figure, then the whole segment joining P and Q must be contained in 7^. Moreover, P 



is a particular polygon whose sides are arcs of circumference. 
'We say that an arc belongs to a circle when it corresponds to a portion of its circumference 



and Q both belong to X^"-~^\ii, i2, k), which is also a convex set (since it is the intersection of convex 
sets), and therefore the segment V=PQcI^^~^^ {ii, 12, k), and also ycX("~^)(ii, i2)- Since PE'ji^ and 
Q^7i2, then VD'ji^^^, and analogously Vn'ji^^fl). It follows, from the fact that ycX'^"~^)(ii, ^2), that 
V cannot hit any circumference other than those of and 742- In addition, the point of intersection 
between V and the circumference of 7^^ must belong to the arc of this circumference contained in 
This arc is ar, since and at are two non consecutive arcs of A. Hence, there is one point 
of Or which belongs to V. As a side of A, arCl^'^~'^\k), whereas Vcyp this point of intersection then 

(n) 

belongs to Xj^^ which therefore is a non-empty set. The same holds for the point of intersection 
between V and the circumference of 7^3 and, consequently, for the fraction of V between the two 
intersection points. 

Finally, we observe that the found set has non-zero measure by construction, and therefore 
0(x|^'' . |)>0. Note that this proof does not hold for n=A, where all the arcs of A are consecutive. 
However, in this case the thesis can be proved in an analogous way. 

In the third case, there is a circle jj whose associated set of arcs a.{j) contains at least two arcs of 
A. It can be shown that two consecutive arcs cannot belong to the same circle. Once observed this, 
the proof is analogous to that of the previous case by choosing at and as non-consecutive arcs 
belonging to Q:(j). ■ 

The previous theorem represents a powerful tool for testing the existence of the intersection regions 
generated by n circles, once those generated by n—1 are known. In fact, the intersection of a set S 
of n>4 circles is non-empty if and only if all the intersections of the subsets of n—1 circles of S are 
non-empty. This simple principle substitutes involved and time-demanding geometric considerations. 

The following theorem states an important property which can be used to derive an effective way to 
compute the area of any intersection of n circles as a function of some specific intersections generated 
by n—1 and n—2 circles. 
Theorem 2 

Consider a subset of n circles 5={7i^, . . . , 7i„}CC, with A<n<Nc, and the associated intersection region 

i \, with </)(xl"^ , i)>0. Denote with m the number of arcs of circumference that delimit 
Then, ifm > 4, there exist two circles 74, 7, 6 5 such that: 

Pi 7i C 7t U 7, (4) 
Proof: See Appendix |Al ■ 



Theorem [2] requires A to have a number of arcs greater than 3. If this does not hold, X*^") collapses 
into the intersection of two or three circles, which can be obtained geometrically. If four or more 
circles are involved. Theorem |2] ensures that the area of X*^") can be found by algebraic manipulation 
of the intersections among n—1 and n—2 circles. 

In fact, if Theorem [2] holds, the intersection area X*^"~^^(t, r) among all circles but 7^ and 7,. is fully 
included in the union between 7^ and 7^. The area of this intersection can be calculated by considering 
the partition of 7iU7r into three regions, one belonging only to 7^, one to both 7^ and 7^, and the 
third only to 7^. All these three regions exist, the second one by hypothesis, the first and the third 
one because otherwise the required intersection X*^"^ would be fully included in one circle, and the 
problem could be reduced to an analogous one with a lower number of circles. Call these areas Bi, 
B2 and Bs. X'^"~^)(t, r) intersects B2 by hypothesis, and it may also intersect Bi, B^, or even both. 
Assume that it intersects all three areas, and call the intersections respectively mathcalAi, A2 and 
^3. It is clear that A2 = These three areas are all unknown. However, it follows from their 
definition that AiUA2= X("-i)(t), ^2 U ^3 = X("-i)(r), and ^1 U A U ^3 = X("-2)(^^ j^en 

0(A) = (0Ml)+</)M2)) + (0M2) + 0(A)) - (0(A)+0(A) + 0(A) (5) 

and, therefore: 

0(XW) =0(x("-i)(t)) + 0(j("-i)(r)) -0(x("-2)(t,r)) (6) 

It can be similarly shown that even if X*^"^^)(t, r) does not intersect Bi, B^ or both, the result still 
holds. 

III. The Algorithm 

In this Section, we describe in detail the proposed algorithm. We remark that the algorithm computes 
both the non exclusive intersection areas (set X) and the exclusive ones (set £). Since the cardinality of 
X grows exponentially with the number of circles, it follows that also the complexity of the algorithm 
is exponential in the worst case (the cardinality of £ instead increases with the square number of the 
deployed circles). Indeed, it is shown in Section |IV] that this happens when a common intersection 
exists among all the circles, with no nested circles. In this case, the overall complexity grows as 
iV^2^'=. The algorithm, exploiting the previously derived results, iteratively checks the existence and 
computes the intersection areas thus achieving a greatly reduced complexity in most cases. 



As stated before, the algorithm is based on a trellis structure, which is built iteration by iteration 
with a simple procedure. After having retrieved all the elements of X, the exclusive intersection areas 
are also derived and collected in a set of vectors, which is the output of the algorithm. 

In this Section, we first describe the trellis structure, showing why it is useful to represent and 
order the intersection areas, and how it is built through simple transition matrices; secondly, we list 
the auxiliary variables which are computed at each step of the algorithm, and explain how they are 
used; finally, we show how the auxiliary variables are updated at each step, and how, at the end of 
the last step, the vectors containing the exclusive intersection areas can be retrieved from them. 

A. Trellis Structure 

The algorithm is based on a trellis structure. In each iteration, the algorithm takes as input the 
structure and the variables built at the previous steps, in order to update the trellis. 

Let us clarify how the algorithm works via a graphical example. Consider the circles in Fig. |3l 
where the first four steps of the algorithm are depicted in the subfigures A, B, C and D, respectively. 
In the first step, the algorithm simply calculates the areas of all the circles. In the second step, the 
algorithm checks the existence and computes the intersection areas of any pair of circles via simple 
geometric calculations. In the third step, the intersection areas of any existing intersection between 
any triplet of circles are computed. This is the last step involving geometric computations, as in the 
subsequent steps the existence and the measure of the intersection regions of any set of 4:<n<Nc 
circles are carried out via simple algebraic manipulation of the areas computed at the previous steps, 
exploiting the two theorems presented in Section |ll] 

In order to effectively keep track of the existence and value of the regions computed in each step we 
build a proper graph, in which the various intersection areas correspond to the vertices of the graph. 
In this way, the relationships can be represented as the edges of the graph. As a byproduct, the graph 
is also useful to list and order all the areas to be computed. 

Consider a set C of Nc circles, labeled as 7,, with l<i<Nc. Our trellis structure is given by a set 
of vertices V and a set of edges TZ. Each vertex corresponds to an intersection among some of the 
circles belonging to the set C considered. Hence, there is a one-to-one relationship among the vertices 
and all possible subsets of C, except the empty set. We define as S{v)CC the set of circles whose 
intersection corresponds to vertex v in the graph. 

The vertices are divided into the subsets V{n), with n=l,2, . . . , N^. The subset V(n) is the set of 



all the intersections among n circles out of the considered (that is, the set of all the subsets of C 
with exactly n elements). 

The obtained subsets can be ordered from V(l) to V{Nc). Each edge can connect only two vertices 
belonging to two consecutive subsets V{n) and V{n+1). An edge connecting jEV{n) and iEV{n+l) 
exists if and only if S{j)cS{i). In other words, i is reachable from j if it corresponds to the intersection 
of all the circles associated with j plus another one. 

Finally, also the elements of each V{n) can be ordered. Vertex i of the graph can be uniquely 
identified by a binary sequence of Nc elements hi=[b\,b2, ■ ■ ■,b^iyj such that bl=l if 'jt^S{i). For 
instance, if Nc=5, the vertex corresponding to the intersection of 72, 73 and 75 can be labeled as 
01101. An equivalent labeling is obtained using the decimal representations of the binary sequences. 
In the example above, the same vertex is hence labeled as 13. The vertices of V{n) are then ordered 
with decreasing labels. An example of the full graph for Nc=5 is also reported in Fig. S] 

The idea behind the trellis is straightforward: a vertex belonging to V(3), for instance, corresponds 
to the intersection of three circles. It follows that this area is related to the three intersection areas 
among any pair of the three considered circles (represented by three vertices belonging to V2), since 
it can be derived from each one of them by adding the missing circle. The exact way through which 
this can be done is explained below, and relies on Theorem |2] 

With the ordering described above, the graph can be fully described by Nc—l transition matrices. We 
define the transition matrix M^^|^, for l<n<Nc, as a binary matrix containing the information about 
which edges exist between the vertices in V{n) and the ones in V{n+1). M|f^|]^(z, j) = 1 if there 
exists an edge connecting the i-th vertex of V{n+1) and the j-th vertex of V{n), and otherwise. 

The computation of the transition matrices can be performed recursively, if we consider an additional 
subset V(0), which contains a single element corresponding to the empty set 0, that is, to the region 
not covered by any circle. This "virtual" vertex is not included in the trellis in Fig. |4l though it can 
be inserted on the left side and connected to all the vertices belonging to V(l). 

With this modification, the following properties about the transition matrices hold: 

• M^^^l^, for 1 < < Nc, is a p x g matrix, where 




(7) 




is a A'^c X 1 column vector whose elements are all equal to 1 . 




where indicates the transposition along the 



secondary diagonal, i.e., the element of the matrix of the i-th row and j-th column is moved to 
the g— j+l-th row and p—i+l-th column. 
• Each transition matrix can be derived recursively as: 



^(^=-1) I 

M^^^"'^ 



(8) 



where I is the identity matrix. We point out that, in any case, it must be n < A^^^. However, 
for n>\Nc/2], due to the third property of the transition matrix, we can compute instead 
^N^^^-n-i Nc-n- ^hc rccursion appears in the first and in the fourth term of the matrix on the 
right side in Eq. ([8]). In the first term both Nc and n are reduced at each iteration, until we get to 
a matrix with n = 0, which can be derived using the second property of the transition matrices. 
As regards the fourth term, at each recursion only A^^^ is reduced. In this case, the recursion ends 
when Nc — 1 = n + 1. \n fact, with the exception of the trivial case n = 0, the matrix M^^;^^^'' 
can now be derived from Mq^""^'', due to the third property, which is in turn computed according 
to the second property. 

Different transition matrices can be multiplied together, thus giving information about which vertices 
of a set V(n) are connected to vertices belonging to V{n+k), with k>l. The elements of the transition 

matrices are binary variables, and we define the transition matrix between V(n) and V{n+k) as 

fc— 1 

B. Auxiliary Variables 

As stated before, the algorithm at each step n finds the intersection areas represented by the vertices 
belonging to V{n). In order to obtain these areas, it exploits the areas belonging to V(n— 1), found 
in the previous step, and applies Theorem |2] based on the relationships expressed by the transition 
matrix M^^\^„^ During the subsequent steps, however, some information must be collected and stored. 
Before describing the steps in detail, we list the auxiliary variables which are computed at each step. 
At the end of the algorithm, they are used to retrieve the required intersection areas. 

• We define the n-th label vector L„ as the vector containing the decimal labels of all the elements 

/ATX 

of the set V(n). Therefore, L„ has length ( j , and its elements are sorted in decreasing order. 



*In the following we drop the superscript of the transition matrices for notation clarity. 



These vectors are computed using the transition matrices, since Li = 
and 

Ln = ^-M„_i,„L„_i (10) 
n — 1 

• The n-th existence vector E„ is a binary vector containing information on the intersection areas 
of n circles; therefore, its size is equal to the cardinality of V{n). The i-th element of E„ is equal 
to 1 if the intersection among the circles associated with the i-th vertex in V{n) is not empty. 
Assuming that all the Nc circles belonging to C have positive radius, then Ei is an NcXl column 
vector whose elements are all equal to 1 . The other vectors can be recursively computed with the 
following rule: 

E„+i = max (M„,„+iE„-n 1, 0) (1 1) 

where 1 and are here column vectors of length I ^ ], with all elements equal to 1 and 



respectively, and the maximum is taken element-wise. Equation (fTTI) can be explained as follows. 
The intersection area Au, of n+1 circles, represented by vertex wEV{n+l), can exist only if all 
the intersection areas among any n of those circles also exist. These areas are in turn represented 
by n+1 vertices belonging to V(n), which are all connected to w in the graph. Assume that w 
is the z-th vertex of V{n+1): the existence of Aw is indicated by En+i{i). Note now that the 
multiplication of the i-th row of M^^+i and E„ is equal to the number of existing intersections 
among n circles out of the n+1 associated with w. Therefore, in order for A^ to exist, this 
product must be equal to n+1. This holds for each element of E„_|_i, which is then reduced to a 
binary vector by subtracting n and nulling the negative elements. When we are considering sets 
of up to two circles we need to do an additional check. In particular. This means that even when 
E„+i(i) = l, the corresponding area may not exist, and the existence vector must therefore be 
modified accordingly. This check is not needed for sets of three or more circles due to Theorem 

m 

• We define the n-th area vector A„ as the vector containing the values of all the intersection areas 
among n circles. As L„ and E„, it is a column vector of I ^ ] elements. The values corresponding 
to non existing areas are equal to 0. We remark that the outcome of these calculations is the set 
X''"\ n=l, . . . , Nc, while our goal is to compute also the exclusive intersection regions, i.e., the 
set £^"'\ We will show how to compute those areas later. 



Besides those vectors, the algorithm also keeps track of the vectors r, and y^, containing the 
radii and the coordinates of the centers of the n circles. Finally, the symmetric matrix D contains the 
distances among the centers: D(z,j) is the distance between the centers of the circles 7^ and 7^. 

C. Computation of the Intersection Areas 

Based on the trellis structure and the auxiliary variables introduced above, we now illustrate how 
the algorithm works. 

We recall that the aim of the algorithm is to calculate the areas of all the exclusive intersections of 
the set £. The algorithm uses the auxiliary variables L„ and E„ to efficiently compute the elements 
of the vectors A„, for n G {1, 2, . . . , A'^c}- 

The algorithm is divided in three phases: 

• Initialization phase, where all the auxiliary variables are initialized; 

• Trellis Computation phase, whose aim is to compute the values of all the vectors A„, for n E 
{l,2,...,iVj; 

• Data Processing phase, where the exclusive intersections of the set £ are retrieved from the vectors 

An. 

We report in the following how the algorithm works in each phase. 

In the initialization phase, the auxiliary variables Li, Ei and Ai are calculated as stated in Section 
IIII-B[ Also, the transition matrices are recursively derived. We assume that all the Nc circles involved 
have finite and positive radius, and thus the measure of the region covered by each circle is strictly 
positive. 

The Trellis Computation phase then is performed in Nc — I subsequent steps. At each step n, the 
aim of the algorithm is to calculate and store the corresponding area vector A„. To this purpose, it 
performs the following operations: 

• generation of L„ from L„_i using (flOl) . This variable is useful, since the binary labels, which can 
be easily obtained from the decimal ones, provide an effective way to recognize which vertices 
of V(n— 1) are connected to each vertex of V(n); 

• generation of E„ from E„_i using (fTT)) . The number and positions of the nonzero elements 
of E„ provide information on which intersection areas must be calculated. For n<3, however, a 
geometric check is necessary for each of these elements, since a nonzero value does not necessarily 
mean that the intersection represented by the corresponding vertex exists. The check is done based 



on the centers and the radii of the involved circles (identified through the label vector L„), and 
it is then possible to obtain E„. No checks are needed for n>3, thanks to Theorem [Tl 
• calculation of A„ from the already known vectors Aj, with i < n. If n < 3, the existing areas, 
according to the information given by E„, can be calculated geometrically, applying known 
formulas with the centers and the radii of the intersecting circles, again retrieved through L„. 
If on the contrary n>3, then we proceed as follows. For the fc-th element of A„, that is, A„(A;), 
we can consider a reduced version of the trellis, containing only the circles whose intersection 
is represented by A„(A;). In this reduced trellis, all the vectors Aj, for i < n, are known, each 
of them containing a suitable subset of the elements of the vectors of the original trellis. The 
problem is now equivalent to the more general one of finding Ajv^ in a trellis when all the other 
vectors Aj are already known, and we refer to this case from now on. 
To solve this problem, we define the vectors Aj, for i G {1,2, Nc — 1}, as: 



Now An^, which is a scalar, can be found from the vectors A;. More specifically, we state that, 
for Nc > 4, Aat^ is equal to the maximum element of —An^^2- This holds also for Nc = 4, 
unless the minimum elements of Ai and A3 are equal, and greater than the maximum element of 
— Ajv^_2- In this special case, an additional geometric check is required, since the value of Ajy^ 
may instead be equal to the minimum element of A^^^i. A formal proof of these statements, as 
well as the required check for the special case Nc = 4, is reported in Section IIII-D[ and is based 
on Theorem |2l 

After the last step of the Trellis Computation phase, all the vectors A,, with i E {1,2, .. . ,Nc}, 
are known. However, the elements of these vectors are not belonging to the set £, since they are the 
areas of the non exclusive intersections. The Data Processing phase of the algorithm performs the 
computation of the exclusive intersection areas, contained in the vectors Aj, with i E {1,2,..., Nc}. 

It can be shown (a sketch of the proof is reported in Section IIII-D|) that the following equation 
holds: 




(12) 




(13) 



where the equality between A^v^ and A^^ is intuitive, since the non exclusive and the exclusive 
intersection among all the Nc circles are necessarily equal. Equation ([TJt is similar to (fT2)) . except 
that now the sum is up to Nc, since here all the vectors Aj are known. The elements of the vectors 
Aj are the required areas, and this concludes the algorithm. 

D. Proof of algorithm correctness 

We prove here that, if the vectors Aj and Aj are known, for 1 < i < Nc, then A^^, which is a 
scalar, can be found as the maximum value of —A]\f^_2- This holds for Nc > 4, and very often also 
for Nc = 4: where, however, in some special cases the required value is instead equal to the minimum 
value of A]y^_i. Finally, for Nc < 3, all the areas can be found geometrically. 

We recall that the vectors Aj contain the non exclusive intersection areas, whereas the vectors Aj 
can be obtained using (fT2l) . The proof is structured as follows: 

• we first determine an expression for the elements of the vectors Aj. We focus on the first element 
of Ai, since all the other ones can be retrieved in an analogous manner. We express them as sums 
of exclusive intersection areas. 

• using the expression of Aj and applying Theorem |2l we prove the statement for Nc > 4. 

• we point out in which cases the statement is not true for Nc = 4, and determine how the correct 
value of Ajv^ can be found in these special cases. 

In the previous Section, it has been stated that the two geometric Theorems [T] and |2] cannot be 
applied directly in the proposed algorithm. This is because when the algorithm is executed, it is 
not known a priori how the circles intersect with each other (which would require an exponentially 
complex conditions check). The only available data at the i-th step are the numerical values of the areas 
computed in the previous i — 1 steps. More precisely, at step Nc, the vectors Aj for 1 < i < Nc are 
available. However, the elements of these vectors are not the exclusive intersection areas, as explained 
before. Vectors Aj, for 1 < i < Nc, can also be computed, according to (fT2]) . but also these vectors do 
not contain the values of the exclusive intersection areas, since A^^ is not known, and consequently 
the sums in (fT2l) are up to Nc — I. 

Therefore, it is not straightforward how the value of Ajv^ can be obtained starting from the vectors 
Aj and Aj, for 1 < i < Nc- We focus on A^v^ since, in order to compute each element of Aj, for 
3 < i < Nc, it is sufficient to run the algorithm while considering only i circles. In this reduced 
version of the algorithm, Aj has only one element, and its role is exactly the same as An^ in the non 



reduced version of the algorithm. Finally, the elements of Aj for i < 3 can be derived via geometric 
computation. 

The unknown value A^^ is the intersection area of all the Nc circles. In order to explain how this 
value can be retrieved, we first clarify which are the areas corresponding to the elements of the known 
vectors. Each element of Aj is the non exclusive intersection area of i circles. In other words, it is the 
value of the area whose points belong to all the i considered circles, but which may be included also 
in other circles. In order to distinguish them, we call /ij^^ij, element of A^ corresponding to the 
non exclusive intersection of circles 7jj, 'ji^..., 7^^.. In Figured we depict an example of deployment 
of four circles. In this scenario, for example, jj,i^2 is given hy g + m + o, whereas /ii,3,4 is equal to 
/ + o. We simply call yU the intersection of all the A^^ circles, which is the only element of A^^, and, 
thus, the value we want to find. Analogously, we call the exclusive intersection of circles 

7jj, 7i2..., With reference to Figure |5l we have yu]; 2 = g, and /^i 3 4 = /. Note that these values are 
not necessarily contained in any of the vectors Aj or A,. 

We now want to express the elements of the vectors Aj, for 1 < i < Ai'c as a function of the exclusive 
intersection areas, and we start with Ai, which has exactly Nc elements. The rationale behind equation 
(fT2l) is the following. Since Ai contains the whole areas of all the circles, in order to find the areas 
covered by exactly one circle we need to subtract the fractions of these areas which are shared with 
other circles. We first subtract the areas shared by two circles, which are contained in A2. In doing 
so, since A2 does not contain the exclusive intersection areas, we are subtracting the areas shared by 
three circles more than once. Hence, we need to sum them again. They are contained in A3. For the 
same reason, we then need to subtract again the areas shared by four circles, contained in A4, and so 
forth. If we had the intersection area of all the circles, at the end of this sum we would get the 
exclusive intersection areas of one circle, that is, with a slight abuse of expression, the areas covered 
by exactly one circle. 

For the sake of clarity, we recall again the example of Figure [5l The first element of Ai contains 



the area of 71, which we call /ii. According to equation (fT2l) . the first element of Ai is computed as: 

Al(l) = /il — /il,2 — ^^1,3 — A*l,4 + /^l,2,3 + /^l,2,4 + /^l,3,4 

= {a + e + f + g + l + m + o)-{g + m + o) + 

-(/ + / + m + o) - (e + / + o) + 

+ (m + o) + (o) + (/ + o) 

= a + o (14) 

As can be observed, the result is not the area covered only by circle 71, which would be instead if 
A4 = o was subtracted by the result. Since A4 is the area we are looking for, it is necessary to find a 
way to determine its value. 

We first derive an expression for the elements of the vectors A^. It is clear from (fT2]) that the 
number of elements of is equal to the number of elements of for each k. Since the derivation 
is analogous for every value of k, we focus on Ai. 

The first element of this vector is obtained as: 

Afc-l 

Ai(l) = /ii+^(-l)"~^ /iMi,i2,...,i„_l 
n=2 l<h<l2<...<in-i<Nc 

= 5^(-ir'A„ (15) 



n=l 



where we define 



Ai — /ii; 



An — ^ fJ'l,il,i2,...,in-l (16) 

l<n<j2<...<i,i_i<Afc 

Since we are computing the first element of Ai, we are focusing on circle 71. The first term of the 
sum is the area of the circle. The second one, A2, is the sum of all the non exclusive intersection areas 
between 71 and another circle, A3 is the sum of all the non exclusive intersections between 71 and two 
other circles, and so on. These areas are not disjoint, since they are the non exclusive intersections. 
We want to rewrite them in terms of disjoint areas, in order to simplify the expression in (fTST ). Define 
the following sums of disjoint areas: 

A* * 
1 = 1^1 

K = Y /^i,ii,i2,...,i„-i 

l<n<j2<...<i„-i<Afc 



According to the above definitions, is the area covered only by 71, is the sum of all the areas 
covered only by 71 and another circle, that is, the sum of all the exclusive intersection areas between 
7i and another circle. With reference to Figure |5] we can write for instance: 

Xi = a + e + f + g + l + m + o , A* = a, 
X2 = e + f + g + 2l + 2m + 3o , \l = e + f + g. 

We can now express the non exclusive intersection areas, which appear in (fTSi) . as functions of the 
exclusive ones, by writing relationships between the Aj's and the A*s. In order to do so, we note the 
following facts: 

• Ai is the area of 71, and is hence given by the sum of all the exclusive intersection areas involving 
7i, that is: 

Ai = 5^ A* (18) 

1=1 

• A2 is the sum of all the non exclusive intersection areas between 71 and another circle. In doing 
this sum, the areas covered by more than two circles (one of which is, by definition, 71), are 
counted more than once. More precisely, the areas covered by n circles are counted n — 1 times, 
and hence: 

Nc 

X, = 5^(^-1) A* (19) 

1=2 

• in general, when computing A^, the areas covered hy n > k circles are counted as many times 
as the number of possible extractions of A; — 1 circles out of n — 1. On the contrary, the areas 
covered hy n < k circles are never counted. Therefore: 

Nc ^ ,N 

n=k ^ ^ 

We can order the expressions of the Aj in an Ai'c — 1 x matrix L, whose z-th row contains the 
terms of the sum defining Aj. The sign of the even rows is then changed, to get: 



L = 
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©A3 •• 
-(DA3 •• 


-[ 1 )^Nc 
(Nc-1\ \ * 
[ 2 j'^A'c 
















With this representation, according to (fTST l. the first element of Ai is simply given by the sum of 
all the elements of L. However, while in (1151) the sum is computed row by row, we now observe that 
it is simpler to consider the columns. The sum of the elements of the first column is clearly equal to 
A^. For 2 < n < Nc — I, the sum s„ of the elements of column n is expressed as: 

n— 1 / ^ 

n — I 



/ ^" = (22) 



k=0 ^ 

The sum of all the columns is hence 0, with the exception of the last one, due to the fact that the sum 
in (fT2l) is up to Ai'c — 1. The calculation of s^, is straightforward: 



Nc-2 



k=0 
Nc-l 



k=0 



-if^Xk (23) 



Putting all together, we have that 

Ai(l) = E 



n=l 

Si + Sn, 



= K + i-if^Xk 

= /it + (-l)^>t,2,3,...,7V. 

= /it + (-!)% (24) 

where /i, as defined above, is the intersection area of all the Nc circles. It is clear that the same 
derivation can be done for all the elements of Ai, such that Ai(A;) = + (—1)^'^//. It can also be 
shown that an analogous procedure can be followed to determine the elements of the other vectors 
An, with 2 < n < N^. In this case, for each element A„(A;) the sums Aj and A* can be defined in the 
same manner, for n < i < Nc, and an (Nc — n — 1) x (Nc — n) matrix can be constructed, finally 
resulting in 

An(A:) = + (-1)'^^""+V (25) 



In the expression above, the indices 11,12, ... ,in depend on the element which is being calculated. 
Reporting again the example in Figure |5l we have: 



Ai 







^1,3 


-/i 




-fi 




-/^ 


^'■2,4 




A''3,4 





A''l,2,4 + A* 
A''2,3,4 + 



(26) 



The vectors above can be computed by the algorithm at each step. We stress again that the area /i, 
as well as all the exclusive intersection areas, are unknown, and cannot in general be retrieved from 
the vectors Aj. This is true for > 3, as stated before, since otherwise all the intersection areas can 
be found geometrically. We skip for now the special case A^c = 4. For N^. > 5, we can now exploit 
Theorem |2l The theorem holds only if the number m of circular arcs that delimit the intersection area 
of all the circles (in this case equal to fi) is greater than or equal to 4. 

If this is not true, then there is a circle 7^ that fully contains the intersection area of all the other 
Nc — l circles. This area also belongs to 'j^, which implies that there is at least one exclusive intersection 
area among Nc — I circles which is empty. Looking at the example for Nc = 4, whose vectors are 
reported in (|26] ). this means that one of the values fJ^* j f. in A3 is zero. As a consequence, in order to 
retrieve it is sufficient in general to take the minimum element of Aat^-i. 

If on the contrary the hypothesis of Theorem |2] holds, then it means that there are two circles 7j and 
7a; that fully contain the intersection area of all the other Nc — 2 circles. Following the same reasoning 
as above, it can be concluded that there exists at least one exclusive intersection area among Nc — 2 
circles which is empty, and the value of fi is the maximum of —Aiy^_2- 

Unfortunately, this is not enough, since it is not known a priori whether the hypothesis of Theorem 
[2I holds or not. If it does not, however, we can use the following argument. If m = p, with I < p < 3, 
there are p circles whose intersection is fully included in any other circle (for p = 1, there is one 
circle which is contained in any other circle). We call V the set of these circles, whereas Q is the set 
of the remaining N^ — p circles. Consider Nc—p — 2 other circles belonging to Q. The intersection of 
these circles and the ones belonging to V is fully contained in the remaining two circles of Q. As a 
consequence, at least one of the exclusive intersection areas among Nc — 2 circles is empty, and again 
fi can be obtained as the maximum of —Aiy^^2, as in the case where the hypothesis of Theorem |2] 
holds. The abovementioned consideration clearly requires that Nc is at least equal to 5, otherwise, if 



p = 3, it is not possible to take — p — 2 circles from Q \j. 

The only case to be studied separately is Nc = 4, which is analyzed in depth in Appendix |Bl 

IV. Algorithm complexity 

The computation of the algorithm complexity is not straightforward, given its strong dependence on 
the specific circles deployment. We therefore proceed to determine an upper bound. We first observe 
that the complexity not only depends on the number of circles considered, but rather on how 
these circles are deployed. For instance, the time needed by the algorithm when the circles are not 
intersecting at all grows as A^"^ (linearly for the calculation of the areas, quadratically to verify the 
absence of any intersection). However, the complexity grows when more circles are intersecting. 

We now make the following observation. The algorithm gives, as a result, both the non exclusive 
and the exclusive intersection areas of Nc circles. Both these results can be useful, depending on the 
considered application. The most demanding task regards the non exclusive intersection areas, whose 
number can rise up to This number is reached when an intersection among all the circles 
exists. If this is not true, only a subset of the areas represented by the vertices of the graph reported in 
Fig. |4] must be computed. Due to Theorem [H this subset is immediately identified after the calculation 
of the non exclusive intersection areas among triplets of circles. In the following we then consider the 
case where the intersection between all the Nc circles exists. 

As regards the exclusive intersection areas, instead, their number is no greater than the number of 
disjoint areas in which the plane is divided when the circles are deployed. It can be observed that this 
number can be at most equal to N^ — Nc + 2, which happens when every circumference intersects all 
the remaining circumferencesc. 

The worst case therefore occurs when there exists an intersection among all the Nc circles and no 
circles are fully included in other ones. We will focus on this case from now on. 

We start with the computation of the non exclusive intersection areas. Having already computed 
the ones created by couples and triplets of circles, we still have to obtain those among 4, 5, . . . , /c 
circles. We observe from equation ([T2] ) that in order to find each element of the vector Aj, with 

^if Nc = 5 and p — 3, the intersection of the circles belonging to V is fully included in the two circles belonging to Q, and the 
same reasoning still holds. 

^This can be easily proved by adding a circle at a time. The first circle creates two disjoint areas. When the i-th circumference is 
added, it can intersect at most all the other i — 1 ones already deployed, each one in at most 2 different points. These 2{i — 1) points 
divide the added circumference into 2(j — 1) arcs. Since all the i—1 already deployed circumferences intersect each other by hypothesis, 
each of these arcs divides an existing intersection area into 2 areas, hence creating a total of 2(i — 1) new areas. It follows that, when 
all the Nc circles have been deployed, the total number of exclusive intersection areas is given by 2 + X]ii2 2(* ~ 1) = -^c — A^c + 2. 



4 < i < A; we need to compute the maximum between a subset of the elements of Aj_2» which is 
in turn obtained as Aj_2 — M^2,i-iAj-i. AH the matrices M^^ are (^'') x (^'') binary matrices, with 
exactly (^'i^*) non-zero elements per row. Therefore, they can be more efficiently replaced by lists of 
indices, each containing (^'') x {^''~^) indices. These lists should be precomputed, possibly in a smart 
manner, exploiting the recursive formulation expressed in ([8]), and the symmetry between M„ „+i and 

^Nc-n-l,Nc-n- 

With these matrices available, we focus on the number of operations needed to find the vector Aj, 
with A < i < k. The vector has (^'') elements. The vector Ai_2 has (/^g) elements. The computation 
of each one, following the expression Aj_2 ^ j-iAj-i, requires Nc — i subtractions, since this is 
the number of non-zero elements in each row of i-i- Summing over all the values of i, the total 
number of subtractions is 

E = o^'^ (2^^ - 2iVc - 2) (27) 



1=2 



Having derived the vector Aj_2, each element of Aj is found as the minimum among a subset of (*) 
elements of Aj_2- Assuming that finding the minimum between n elements requires n comparisons. 



the overall number of comparisons needed is: 

E (t) (2) = \ (^""^ - + ^) " ^^^^ 

We have shown that, in order to compute the 2^" non exclusive intersection areas of a family of 
Nc circles, the number of subtractions needed grows as N^fl^", while the number of comparisons 
needed grows as N^2^''. However, in a smart implementation, once an element of Aj is found, it may 
also be compared with the i elements of Aj„i which correspond to the intersections of all but 1 of 

(i) (i— 1) 

the i circles. In fact, whenever X). . . , is equal to X). . it follows that the circle 7j fully 

contains the intersection of the other i — 1 circles, and therefore all the exclusive intersections between 
subsets of circles containing 74^,7^2, . . .ji-^^ and not containing 7j- can be immediately set to zero. 
In this manner, at the end of the computation of the non exclusive intersection areas, several of the 
non existing exclusive intersection areas have been also identified. 

Having computed all the non exclusive intersection areas, the last step is to compute the exclusive 
ones. Instead of using equation (fT3l) , the following relationship is also valid: 

j=i+i ^ ^ 

. Ai i=N, 



which can be used recursively, from A^v^ back to Ai. Although it seems that the same number of 
operations is involved, this is not actually true. This is due to the fact that the number of exclusive 
intersection areas is much lower, as said above. Starting from the exclusive intersection of all the 
Nc circles, which is known (being equal to the non exclusive one), the intersections among A^^^ — 1 
circles are retrieved by means of matrix _^ through the equation defined above. However, some 
of these areas may result equal to zero. This information can be immediately exploited, by deleting 
the corresponding elements in the matrix _2 7v^_i (which can be written in the form of a list, as 
explained above). The result is that when computing the intersection areas among Nc — 2 circles, only 
the required operations are performed, neglecting the subtractions of non existing areas. 

Now, in order to compute the number of required operations, we refer to the symmetric case where 
all the Nc circles have the same radius, and when the intersection between all of them is bounded 
by a circular polygon with Nc sides, each belonging to a different circle. It is easy to verify that by 
changing the positions or the radii of the circles, the number of exclusive intersection areas cannot 
increase any more^. 

In the configuration taken into account, it can be proved that each non exclusive intersection between 
k circles contains (^"^ 2^^^) + 1 exclusive intersection areas. The corresponding exclusive intersection 
area, therefore, can be computed through {'^"^2^^) subtractions. The overall number of operations, in 
order to find all the exclusive intersection areas, is given by: 

g {^k) (^^ ~2 + 3) - ANciNc + 1)) (30) 

which grows as fast as A^'^2^^ 

We conclude that, in the worst case, where all the Nc circles intersect each other and there are no 
nested circles, the complexity of the algorithm, in order to find all the 2^" non exclusive intersection 
areas and all the A^"^^ — Nc + 2 exclusive intersection ones, grows as 2^''Nl. 

A. Comparison with Monte-Carlo approximation 

Other numerical methods can be used to approximate the values of the intersection areas. One of 
the most common ones is the Monte-Carlo approximation: a number Np of points is randomly chosen 
in the area where the circles are deployed. For each of them, the distances from the Nc centers are 

'This number actually does not change, as long as the hypothesis about the intersection area among all the Nc circles holds, except 
for cases where three or more circles intersect in a single point. However, these cases lead to a lower number of exclusive intersection 
areas; in addition, if the circles are randomly deployed, the probability of these configurations is 0. 



computed, thus finding the intersection area it belongs to. Since the probability of choosing a point 
in a given intersection area is equal to the ratio between this area and the total deployment area, 
the values of the intersection areas can be approximated using a high enough number of test points. 
Although the Monte-Carlo method is simpler to implement, we notice that Nc distance computations 
and comparisons are needed for each deployed point. In addition, the precision of this method depends 
on the number of deployed points, especially if the ratio between the average circle radius and the 
side of the deployment area is small. If this is the case, even checking the existence of the intersection 
areas between small circles may require a huge number of points. Indeed, for an area which is p times 
smaller than the total deployment area, in order to find a value with a relative error equal to e with 
probability A, the number Np of points needed can be approximated as follows. 

When Np points are tested, each of them is within the desired area with probability p. Therefore, 
the total number of points which falls within the desired area is a binomial random variable, with 
mean pNp. If Np is high enough, which is true in our case, then the binomial random variable can be 
well approximated with a Gaussian random variable, with mean pNp and variance = pNp{l — p). 
Call this variable X. The condition is then: 



P [(1 - e)pNp < X < (1 + e)pNp] = A 



(31) 



By introducing Y = X/gx —pNp, which is therefore a Gaussian random variable with zero mean and 
unit variance, we can rewrite the condition as: 

^P^p / / epiVp 



P 



<Y < 



y/pNp{l-p) ^pNp{l-p) 

from which, given the simmetry of the Gaussian pdf, we get: 



= A 



Q 



epNp 



1 - A 



(32) 



(33) 



^pNp{l-p) 

where Q{-) is the Gaussian Complementary Cumulative Distribution Function. Solving for Np, we 
finally obtain the result: 



Nr, 



1 — P 

e^p 



1 - A 



(34) 



For p — 0.1, € — 0.01 and A = 0.9, the required points are almost 250000. The proposed algorithm, 
on the contrary, gives the exact values of all the non exclusive and exclusive intersection areas, 
even when the deployment area is much larger than the average circle area (in this case, indeed, the 
probability that all the circles are intersecting with each other is quite small, which further reduces 
the algorithm computational burden). 



V. Network Design Application 



The algorithm presented in this paper may become a useful tool to determine distribution functions 
which would be hard to derive analytically. In this Section we show how this can help in network 
design problems, in order to determine the optimal setting of specific parameters. 

Consider a wireless network where fixed Access Points are distributed in a given area, and a mobile 
terminal whose position is randomly chosen in the same area, with uniform distribution. Assume that 
the transmitted power of the mobile terminal is fixed and equal to Pm, and that a target SNR T is 
required for decoding. With a commonly used approximation, consider a circular area around each 
access point as its coverage area. The radius of this area can be chosen based on the average SNR 
or on the outage probability. Once the channel model is determined, the radius depends only on the 
transmission power Pm- 

It is clear that, if the power is high enough, there are regions covered by two or more access points. 
If we assume some sort of cooperation among the access points, users that are located in these areas 
can take advantage of the spatial diversity by transmitting to multiple access points. The network 
topology has a strong impact on the overall performance: if the considered area is fixed, increasing 
the density of access points makes it easier for a transmission to be decoded by several receivers, but 
also implies an increased network deployment cost. 

In the following example, we assume flat Rayleigh fading, such that, if W is the noise power, the 
SNR at distance d from the receiver is given by the well known equation: 



where A is a pathloss factor, a is the pathloss exponent, and h is the channel fading gain, distributed 
as a complex Gaussian random variable with zero mean and unit variance. The attenuation factor a, 
due to shadowing effects, is here considered constant. The radius R of the coverage circles is defined 
as the distance at which the average SNR at the access point is equal to a given value A. Therefore, 



Alternatively, the radius can be determined based on a given outage probability, with no substantial 
difference in the results. Note also that, since fixed transmit power is assumed, the same analysis could 
be done for the downlink as well, with Base Stations cooperating in the transmission phase. Usually, 
however, the coverage bottleneck lies in the uplink, due to the reduced available power at the mobile 
terminals, thus we focus our attention on this scenario. 



SNR{d) 




(35) 




(36) 



In this Section, we set the parameters as in Table HI To analyze the performance, we compute the 
outage probability in the absence of interference, defined as the probability that the SNR is below the 
target value T. If the node is in the coverage area of two or more access points, we assume that a 
Maximum Ratio Combining strategy is applied. A very robust wired channel among the access points 
is assumed, and the fading is considered independent among different channels. Therefore, the SNR 
of the transmission from a source to n access points can be simply written as 

P (T " 

i=l 

which depends on the distances from all the access points within coverage range. The circumferences 
around the positions of the access points partition the deployment area Ad into J\f disjoint areas. 
Call them Ai, with i E {1, 2, . . . ,J\f}, and call Lp(i) the number of coverage circles Ai belongs to. 
If we call V[Ai] the probability that the source node is located in Ai, the global outage probability 
^ = V[SNR < r] has the general expression: 

1=1 "'O "^0 

(38) 

Note that (^(i) may also be equal to 0, if the chosen parameters and topology do not guarantee 
full coverage of the deployment area. For users located in these areas, we consider that the outage 
event has probability 1 . In general, J\f depends on the selected topology, and grows quadratically with 
the number of circles. We point out that, if interference is also taken into account, the SNR should 
be replaced with the Signal to Interference-plus-Noise Ratio (SINR). In this case, however, also the 
interference term should be averaged, by integrating over the position of the interferer(s). Although 
a simplified interference model may be used, we notice that this model should take into account the 
correlation between interference levels at all the Base Stations connected with the source of the useful 
signal. Such a model is beyond the scope of this paper. 

Equation (|38l) requires the knowledge of three distributions: 

• distribution of the source node position, which is uniform by assumption. This means that V{Ai) 
is given by the ratio between the area of Ai and the area of the whole deployment area Ad- 

• distribution of the SNR, given the number of access points within transmission range and their 
distances from the source node. Once the distances are fixed, it follows from (l37l) that its 
distribution is the distribution of a finite sum of independent exponential random variables, whose 
parameters are related to the distances of the access points. This distribution is known. 



• distribution d ^ ^ ("^i) ^2, • • • , '^i^(j)) of the distances between the source node and the access 

points within coverage range, dependent on the considered area Ai. This distribution strongly 
depends on the shape of Ai. Note that also for areas covered by only one access point, this 
distribution is no longer the distribution of the distance of a point randomly placed in a circular 
area from its center. 

As to the topology, in this example we consider a cellular-like deployment of the access points, 
which are distributed in a hexagonal grid with side equal to L. Every other possible deployment is 
admissible; however, our choice leads to a lower computational burden due to the symmetries of the 
selected topology. In Figure |6l an example of the considered topology is reported, for a fixed L and 
R. The overall performance can be studied as a function of these two parameters, or equivalently, of L 
and Pa/. For the sake of power saving, lower values of Pm are preferable, whereas, in order to reduce 
the network deployment cost, a higher L is desirable. This in turn increases the outage probability, 
thus highlighting the need for a tradeoff. 

The proposed algorithm can be used to find V[Ai] for each area Ai. Moreover, it can also be used 
to find ^2, • • • , as we describe in the following. 

Consider an area Ai covered by only one access point, for example the quasi hexagonal shaded area 
around access point 6 in Figure |6l An analytical expression for the cumulative distribution function 
of the distance from the access point here is hardly derivable. However, we can approximate it by 
running the proposed algorithm several times, and properly selecting the radii of the circles. More 
specifically, we quantize the cdf with arbitarily small step p. To obtain the desired cdf, we keep the 
center of the circles and all the radii, except the radius of the circle around the access point 1, which 
is set to 0, and then increased by p each time the algorithm is run. It is then sufficient to compute at 
each iteration the value of the area covered only by access point 1, and normalize it with the area of 
Ai to get the required cdf and, by numerical differentiation, the corresponding pdf. For areas covered 
by two or more circles, the same method is used, now properly varying the radii of the circles around 
the involved access points. Two and three-dimensional distributions are obtained for areas covered 
by two and three circles respectively. However, for areas covered by four or more circles, only three 
dimensional distributions are to be computed, since three distances uniquely determine the position 
of the source node, and hence all the other distances as well. Therefore, at most three integrals are 
to be numerically computed in (|38T ). Results become more and more precise for smaller values of the 
quantization step p, but the computational burden is also increased. 



If we consider the cellular like distribution of Figure |6l we can assume an infinite deployment area 
Ad- In fact, due to the tessellation of the plane, the analysis can be limited to a finite region (the 
coloured region in Figure |6), whose area can be regarded as the whole deployment area. With this 
choice, border effects are avoided, and a small number of circles are to be considered. The number 
of areas to be computed depends on the ratio between L and R. 

In Figure 111 we report the success probability, that is, V[SNR > F], as a function of the transmission 
power Pm and the distance L between the access points, for F = 10 dB. Once the target SNR F is set, 
and the corresponding graph is plotted, it is possible to determine the required transmission power to 
achieve a given outage probability with a fixed access point density; vice versa, the minimum density 
to achieve the same outage probability can be found when the transmission power is instead fixed. 
Clearly, the success probability is higher when more power is available at the source, or when the 
density of the access points is increased. However, this comes at the cost of a higher number of 
required access points or of a shorter battery life of the source node. Depending on the relevance of 
these two costs, an objective function g{.) could be also defined. Several choices are possible; as an 
example, the objective function may have the following form: 

9{Pm,L) = — — (39) 
uPm + r]V{L) 



where the success probability appears at the numerator, and is derivable through (|38l ). At the 
denominator, 'D{L) is the access points density, equal to 2/(-\/3L^), whereas /i and r] are normalizing 
constants, whose values can be properly selected. As an example, the value of (7 as a function of Pm 
and L is depicted in Figure [8] for specific values of ji and rj. Although the success probability increases 
with Pm and with V{L), the additional cost in terms of energy and number of deployed access points 
penalizes the choice of high Pm and low L, thus identifying an optimal region. 

VI. Conclusion 

In this paper, a practical algorithm for the computation of the intersection areas among any number 
of circles has been presented. The algorithm, based on two geometrical results, is designed to operate 
in an iterative manner, and takes advantage of a trellis structure to order and calculate all the required 
areas, given the radii and the mutual positions of the circles. An application of the algorithm has been 
presented in a network design problem, where cooperation is available among several access points. 
Our algorithm makes it possible to derive the distribution of the location of the source terminal, thus 



allowing the calculation of the outage probability as a function of the transmission power and the 
density of access points. 

Appendix A 
Proof of Theorem [2] 

Proof: Analogously to the proof of Theorem [U we distinguish three cases: 

1) 3{7„j, . . . , 7„^}C5:q;(mi), . . . , q;(m^)=0; 

2) V7,g5 \cx{{)\ = l- 

3) V7,g5 q:(z)^0 and 3j : |a(j)|>l. 

We skip for now the first case, and focus on the following ones. 

In the second case, each circle 7^ has one arc belonging to A. Since m>4, there exist two circles 74, 
7rG5 such that at^(y.{t) and ar^(y.{r) are non-consecutive sides of A. Let us denote as P and Q the 
points of intersection of 7^ and 7^. Note that the existence of P and Q is guaranteed by hypothesis, 
as X(")7^0. We also define as the arc of A belonging to 7^ and fully contained in 7^, and a* as the 
arc of A belonging to 7r and fully contained in 7^. 

Moreover, there exists a circle 'yn^S such that P^'jh- In fact, if this circle did not exist, then PgX*^") 
and the arcs and at would be consecutive arcs of A. 

We need to distinguish two configurations for point Q: 

• In the first configuration, Q does not belong to 'jh (see Fig. |9l). Therefore, the intersection of 
the circumferences of 7^ and 7/j belongs to the arc a* and the intersection of the circumferences 
of 7t and 7/1 belongs to «[. As a consequence, 7/1 is divided into three disjoint areas ^iC7r, 
^2^7^07^ and AsC'jt. We thus conclude that 7/iC7r,U7f. Since, X("~^)(r, t) is contained in 7/j, 
then it is also included in ^r^'jf 

• In the second configuration, Q belongs to 7/1 (see Fig. (TO]). In this case there exists a such 
that Q^jx- If P^lx, the theorem can be proved as in the previous configuration. If PE'Jx, then 
a* and a[ belong to 'jh^lx, as they are sides of A and hence are part of X^"). Thus, the points 
of intersection between the circumferences of 7/1 and 'jx cannot belong to 7rn7(. It follows that 
7^074 is divided into three regions: AiC'jh, A2C'yhf^'yx and A^C'jx- Therefore, 7rn7(C7/iU7i:, 
and since X("~^)(/;,, x)C7,.n7t the thesis is proved. 

In the third case, consider the circle 7^ with |Q;(j)|>l. In the proof of Theorem [T] we showed 
that the two arcs A adjacent to atEa.(j) must belong to the circumference of two different circles, 
say 'jh and 7^. By construction, the arc at is contained in 'jh^lx- Moreover, also arEa(j), with 



ttrT^a*, is contained in ''jh^lx, as is an arc of A. Therefore, 7^ can be divided into three disjoint 
regions analogous to those of the previous case, namely Ai'Z'^x, A2<Z'^x^1h and A^C'^h, which implies 
7jC72:U7/i. Since X*^"^^)(2;, h)C'jj, the thesis in the third case is proved. 

We now conclude the proof with the first case. Here, there exist w circles containing xl"^ . , , with 
l<w<n—3. Define the subset S^^~'^^=S\{'yui, ■ ■ ■ 5 7«„_™}c5 that contains the circles 7^ such that 
a(j)7^0, and their intersection X^'^~'^\ui, . . . ,Un-w)- Thus, there are n—w circles with a non-empty 
set of arcs a{j) and w circles that fully contain I^^^ . |. This case is equivalent to the second or 
third case, if we consider only the circles with a non-empty set of arcs of A. Since we have already 
shown that the theorem holds in those cases, we have here that 37*, 7,.g5"-"' : I^''-'"''^\t,r)c-ftUjr- 
Therefore, 

X("-2) (r, t) C li--'-^) (r, t) C 7.U7t, (40) 

In fact, I^"'~'^\r,t) is equal to the intersection between X'^"^"'~^) (r, t) and nie5\5("-™) ^^^^ 
a subset of x^^^'^^'^\r, t). That proves the theorem in this case. Note that w<ri— 3, in fact in order to 
have m>4 arcs, we need A to be bounded by at least three circles. ■ 

Appendix B 
Derivation of A Nc foR A^c = 4 

We describe here how A can be computed when Nc = A and when an additional geometric check 
is necessary. The available vectors in this case are Ai, A2 and A3, whose expressions have been 
reported in (|26] ). It is useful to consider A2 = — A2. In this manner, all the elements of A2 are smaller 
than or equal to /x, whereas all the elements of Ai and A3 are greater than or equal to /i. Therefore, 
we define 

= min {^i^ ? ^7 = max {^'2^ ? c-y = min {^i^ (41) 

According to the considerations reported above, if m > 4 at least one of the exclusive intersection 
areas between 2 circles is 0, and fi is simply equal to h-y. The same holds also if m = 1 and m = 2, so 
the only case that must be studied is when m is equal to 3. Assume, without loss of generality, that 
the circle 74 contains the intersection of 71, 72 and 73. Then, there are two possible cases, as reported 
also in Figure [TT] 

• 74 is contained in the union of the other three circles. In this case, = 0, and hence /i = a^; 

• 74 is not contained in the union of the other three circles. In this case, it can be shown that it 
fully contains the exclusive intersection of two circles, meaning that /i = 6^. 



Moreover, since the intersection area of 71, 72 and 73 is included in 74, also /Ui 2 3 = 0' both 
cases also jj = c^. In this manner, we have proved that in any case the unknown value of n is equal 
to one of the three values summarized also in Table im However, it is still to be 

determined how the algorithm can recognize which one of the three terms is the actual value. In most 
cases, this can be inferred by the relationships among their values. We can distinguish the following 
exhaustive possibilities: 

• 7^ 6^ 7^ c^: in this case, looking at Table UIl it follows that m 7^ 3, and therefore ji = h^; 

• = ^ or ^ = in both these cases, recalling that a^, > /i whereas < /i, it 
follows that necessarily /i = 6^; 

• = 6^ = c^: as in the previous case; 

• = ^ b^: this is the only case where it is not possible to determine whether jj, = ox jj, = a^. 
In fact, looking at Table UIl this may happen both if m 7^ 3 and if m = 3 and 74 C Ui=i 7*- 
the former case, /i = 6^, in the latter instead /i = a^. 

It is clear that, in the last case, it is very unlikely that m 7^ 3, since this would mean that the smallest 
exclusive intersection among 1 circle has the same (nonzero) area as the smallest exclusive intersection 
of three circles (since = c-y). Anyway, since this may happen, it is necessary to determine the value 
of m in a different way. 

We first of all calculate which values of m may effectively result in = ^ b^. Since it can be 
shown that 1 < m < 2(Nc — 1), for Nc = 4: we have 1 < m < 6. We exclude the following cases: 

• m = 1: in this case 3i: 7j C 'jk, VA; G {1,2,3,4}, k ^ i. This implies that /i* = 0. Moreover, 
Vj 7^ i, we have % D^jj C 7^, V/c G {1, 2, 3, 4}, k 7^ Hence, also = 0, which implies that 

• m = 2: in this case 3i, j, i 7^ j: 7j D'jj C 7^, VA; G {1, 2, 3, 4}, k 7^ This means that = 0. 
In addition, \fk 7^ it is also true that 7,; fl 7^ fl 7^ C E {1, 2, 3, 4}, p 7^ k, which in 
turn implies that also At*j^ = 0. As a consequence, 6^ = c^. 

• m = 4, with 2 sides belonging to the same circle: in this case Theorem |2] holds, meaning that 
there exists one exclusive intersection among two circles with zero area. Moreover, it is clear 
that there exists a circle 7^ that fully contains the intersection of the other three circles (the ones 
whose arcs delimit the intersection area of all the four circles), meaning that u*u„ = 0, with 
j, k,p E {1, 2, 3, 4}, and i ^ j ^ k ^ p. Therefore = c^. 

• m = 5: here again Theorem |2] holds, meaning that 3i,j G {1,2,3,4}, i 7^ j: = 0. In 



addition, as reported in Figure [I2l it is clear that two non consecutive arcs ai and al delimiting 
the intersection area of the four circles belong to the same circle, say 71. Assume that the arc 
between them, namely 0:2, belongs to 72. This arc divides the circle 71 in two parts. The one 
containing the intersection of the four circles is fully contained in 72, since the two circumferences 
7i and 72 cannot intersect in more than two points, and they already intersect in P and Q. The 
other part of 71 is contained in 73 (as well as in 74). In fact, if this were not true, since P and Q 
both belong to 73, the arc (3i should intersect the circumference 73 in two points. Since 71 and 
73 already intersect in R, this would cause them to intersect in more than two points,which is 
not possible. Therefore, 71 C 72 U 73, and = 0, meaning that = b^. 
• m = 6: the same reasoning as for m = 5 can be done in this case, again resulting in = h^. 
Having excluded all the above listed cases, there are only two possible deployments that may result 

in = h-y-. either m = 3 or m = 4 with all the four arcs belonging to different circles. They are 

reported in Figure [13] 

Recall that if m = 3, the only possible case is the one reported in Figure [TT] on the left (and 
depicted also in Figure [l3]on the left). The deployment depicted in Figure [TTIon the right is instead 
not possible, since it implies reported also in Table III 

From the previous investigation, we can conclude that when = 7^ h^, the value of /i is if 
the circles are deployed as in Figure [13] on the left, and is equal to if the circles are deployed as 
in Figure [13] on the right. No other deployments are compatible with the given inequality among a^, 

and c^. The straightforward way to distinguish between the two cases is to calculate all the twelve 
intersection points between the four circles. In both situations each circle contains exactly three points 
of intersection between the other three circles. More precisely: Vz G {1, 2, 3, 4}, 3P, Q,R: P e ^jCi^k, 
Q G 7j n 7p and G 7^ fl %, with j,k,p e {1,2, 3, 4}, and i, j, k and p all different from each other, 
such that P,Q,R G 7^, where again 7^ indicates the circumference of circle 7^. 

However, if m = 3 there is one circle that contains the intersection area of the other three circles, 
meaning that there is one (and only one) circle 7^ such that, using the same notation as above, P G 7p, 
Q G 'jk and R G 'jj. This is not true when m = 4, since in that case, one of the three points contained 
in each circle 7^, given by the intersection of two circumferences ■yj and 7^, does not belong to the 
third circle 7^. Simple geometric comparisons among the distances between centers and intersection 
points and the radii of the circles are then enough to distinguish the two cases and, finally, determine 
the correct value of /x. Note that this check is necessary only for Nc = 4, and only when = ^ b^. 
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Fig. 2. Intersection of n circles. Here each circumference contains one arc belonging to the circular polygon A, which delimits the 
intersection. 



Fig. 3. Example of circles configuration and areas considered in the first four steps of the algorithm. 





Fig. 6. Example of plane tessellation for \/3/3-£/ < R < V3/2L. The shaded area around Access Point 3 represents the basic element 
of the tessellation. The shaded area around Access Point 6 is instead the portion of the area of each cell which is covered by only one 
Access Point. 
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Fig. 7. Success probability, as a function of Pm and L. Here F = 10 dB. 
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Fig. 8. Values of the objective function g, as reported in i39\ . as a function of Pm and L. Here F = 10 dB, fi — 1W~^, and 
T] = 30 




Fig. 9. First configuration of ttie two possible when m=n, and each of the n intersecting circles has exactly one arc of circumference 
belonging to the polygon A. In the figure A2 contains A. 



Fig. 10. Second configuration of the two possible when m=n, and each of the n intersecting circles has exactly one arc of circumference 
belonging to the polygon A. 




Fig. 12. Intersection of 4 circles, with m = 5. 



Fig. 13. The only two possible deployments of 4 circles such that = c-^ =^h-y. 



TABLE I 

System parameters for the cellular-like topology. 



Noise power W 


-lOSdBm 


Path-loss exponent a 


3 


Shadowing margin a 


10 dB 


Fixed attenuation parameter A 


30 dB 


SNR threshold A for decoding radius 


10 dB 


SNR threshold T for outage probability 


10 dB 


Quantization step p 


R/50 m 



TABLE II 

Relationships among a-y, b-y and Oy, depending on the positions of the circles. 



Case 


Ordering 


Value of n 


m 3 




bj 


m = 3 and T4 C U^^-i 7' 


a-, = C-, > h-, 


<i-X= c,) 


rn = :') and - 4 ^ Ui-i "1 


(/,-,, > b~ = (.:-,, 





